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Abstract. I discuss quarkonium suppression in equilibrated strongly interacting matter. After a 
brief review of basic features of quarkonium production I discuss the application of recent lattice 
data on the heavy quark potential to the problem of quarkonium dissociation as well as the problem 
of direct lattice determination of quarkonium properties in finite temperature lattice QCD. 
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1. Introduction 

The behavior of the heavy quarkonium states in hot strongly interacting matter was pro- 
posed as a test of its confining nature, since a sufficiently hot deconfined medium will 
dissolve any binding between the quark-antiquark pair [1]. Another possibility of disso- 
ciation of certain quarkonium states (subthreshold states at T = 0) is the decay into open 
charm (beauty) mesons due to in-medium modification of both the quarkonia and heavy- 
light meson masses [2,3]. 

The production of J/t/j and T mesons in hadronic reactions occurs in part through pro- 
duction of higher excited cc (or bb) states and their decay into quarkonia ground state. 
Since the lifetime of different subthreshold quarkonium states is much larger than the typ- 
ical life-time of the medium which may be produced in nucleus-nucleus collisions their 
decay occurs almost completely outside the produced medium. This means that the pro- 
duced medium can be probed not only by the ground state quarkonium but also by different 
excited quarkonium states. Since different quarkonium states have different sizes ( binding 
energies ), one expects that higher excited states will dissolve at lower temperature than 
the smaller and more tightly bound ground states. These facts may lead to a sequential 
suppression pattern in J/ ip and T yield in nucleus-nucleus collision as the function of the 
energy density [4] . 

The rest of the paper is organized as follows. In section 2 I briefly review the basic 
features of quarkonium production. In section 3 I discuss the possibility of quarkonium 
dissociation below deconfinement. In section 4 I use the lattice data on the heavy quark 
potential to determine the dissociation temperatures due to color screening. Finally the 
problem of direct lattice determination of quarkonium properties at finite temperature is 
discussed in section 5 followed by the conclusions presented in section 6. 
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2. Quarkonium production and feed-down 

It is well known that J/ip production in hadron-hadron collision is due to a considerable 
extent to the production and subsequent decay of higher excited cc states [5-7]. The feed- 
down from higher excited states was systematically studied in proton-nucleon and pion- 
nucleon interactions with SOOGeV incident proton (pion) beams [7]. In these studies the 
cross sections for direct production of different charmonium states (excluding feed-down) 
were measured. Then making use of the known branching ratios B[xi(lP) — > ip(lS)] = 
0.27 ± 0.02, B[ X 2{IP) -> i>0-S)] = 0.14 ±.01, and B[^(25) -> ip(lS)] = 0.55 ±0.05, 
one obtains the fractional feed-down contributions /; of the different charmonium states to 
the observed J/ip production; these are shown in the second and third columns of Tab. 1 



state 


fi(n-N) [%] 


fi(pN) [%] 


state 


Mpp) [%] 


if RQCD (PP) [%] 


JMis) 


57 ±3 


62 ±4 


T(15) 


52 ±9 


52 ± 34 




20 ±5 


16 ±4 


Xb(lP) 


26 ±7 


24 ±8 




15 ±4 


14 ±4 


T(25) 


10 ±3 


8±7 


iP(2S) 


8±2 


8±2 


Xb(2P) 


10 ±7 


14 ±4 








T(35) 


2 ±0.5 


2±2 



Table 1: Feed-down fractions for from higher excited states to the J/ip and T states. 



In the case of bottomonium the experiment provides only the inclusive (i.e. including 
also the feed-down from higher states) cross section for different (nS) states [8]. The 
feed-down from (nP) states is known only for transverse momenta pt > 8GeV/c [8]. To 
analyze the complete feed-down pattern, we thus have to find a way to extrapolate these 
data to pt = as well as to determine the direct cross section for different (nS) states. This 
can be done using the most simple and general model for quarkonium production, the color 
evaporation model [9]. In particular this model predicts that the ratios of cross sections 
for production of different quarkonium states are energy independent. This prediction 
was verified for a considerable range of energies [10]. The ratios between the different 
Xj(1P) states in this model are predicted to be governed essentially by the orbital angular 
momentum degeneracy [9]; we thus expect for the corresponding cross-sections 

X o(lP):Xi(TP):X2(lP) = l:3:5. (1) 
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From Table 1 wehavefor7r~A r collisionsx2(l-P)/xi(l-P) — 1.44±0. 38 and thus reason- 
able agreement with the predicted ratio 1.67. Actually, for pN interactions, the experiment 
measures only the combined effect of xi and \2 decay (30 % of the overall J/tp produc- 
tion); the listed values in Tab. 1 are obtained by distributing this in the ratio 3:5. 

Using considerations based on color evaporation model, in particular Eq. (1), the feed- 
down from higher excited bb states to T production can be predicted [11]; the feed-down 
fractions are summarized in Tab. 1 . Alternatively the feed-down fraction from higher ex- 
cited bb states can be predicted using NRQCD factorization formula [11]. The results of 
this analysis are summarized in the last column of Tab. 1 

3. Quarkonium dissociation below deconflnement 

Recent lattice calculations of the heavy quark potential show evidence for the string break- 
ing at finite temperature [12]. On the lattice the potential is calculated from the Polyakov 
loop correlator, to which it is related by 



where L(r) is the Polyakov loop (see e.g. Ref. [12] for definition). The normalization 
constant C contains both the cut-off dependent self-energy and the entropy contributions 
-TS (for T ^ 0, - In < L(r)L t (0) > is actually the free energy of the static QQ pair). 
For a properly chosen normalization constant C, V(r, T) is the ground state energy of an 
infinitely heavy QQ pair. In absence of dynamical quarks (quenched QCD) V(r, T) is 
linearly rising with r for large separations indicating the existence of a flux tube (string). 
If dynamical quarks are present the flux tube can decay (the string can break) by creating 
a pair of light quarks qq from the vacuum once V(r, T) is larger then twice the binding 
energy of a heavy-light Qq (qQ) meson 1 . Thus the potential at very large distances is 
constant (T) and is equal to twice the binding energy of a heavy-light meson. 

It is expected that medium effects are not important at very short distances. Therefore at 
very short distances the potential V(r, T) should be given by the Cornell potential [14] 



We use this fact to determine the normalization constant C and set the potential to be of the 
Cornell form at the the smallest distance rT = 0.25 available in lattice studies of [12], with 
e = 0.4 as expected for (2+l)-flavor QCD [15]. The resulting potential and Voo(T') are 
shown in Fig. 1. Note the strong temperature dependence of Voo(T). Since for sufficiently 
heavy quarks (toq 3> Aqcd) it does not matter whether the quark is infinitely heavy 
or just merely heavy, the open charm (beauty) meson masses are approximately given by 
2M D . B {T) = 2m c . b + V 00 {T). 

Now the temperature dependence of the different quarkonium states should be ad- 
dressed. At zero temperature the heavy quark masses permit the application of potential 



1 Similar phenomenon occurs of course at T — 0. However it is much more difficult to observe it 
on lattice (see e.g. [13]). 



V(r, T) = - In < i(r)L t (0) > +C, 



(2) 



V(r) = — - + err 



(3) 



r 
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Figure 1. The heavy quark potential and its asymptotic value below deconfinement at 
different temperatures. The line on the bottom figure is a fit to the data points. 



theory for the description of quarkonium spectroscopy (see e.g. [13]). Furthermore it turns 
out that the time scale of gluodynamics relevant for quarkonium spectroscopy is smaller 
than (m,Qv)~ (v being the heavy quark velocity) [13]. For sufficiently heavy quarks this 
time scale is much smaller than the typical hadronic time scale Aq C£) <~ 1/m. The decay 
of the flux-tube like all other hadronic decays has time scale of order 1/m. Therefore 
in the potential theory the potential must always be of Cornell form (i.e. linearly rising 
at large distances). These considerations have direct phenomenological support. Namely, 
simple potential models with linearly rising potential can describe reasonably well also the 
quarkonium states above the open charm (beauty) threshold. Many of these higher excited 
states have an effective radius of order of or even larger than 1/m [15,14]. Contrary to 
this situation in the case of the potential becoming flat around 1/m (the expected radius of 
string breaking at T = 0) the higher excited states above the open charm (beauty) threshold 
simply do not exist. Therefore the temperature dependent quarkonium masses have been 
determined from the Schrodinger equation 



2m Q - —V 2 + V(r) 
m Q 



$i = M i <S> i , (4) 
where the potential V(r) is identified with the temperature dependent string potential [16] 



4 



Pramana - J. Phys., Vol. X, No. Y, March 2002 



3.8 




2.6 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

T/T c 





10.6 






10.4 


3S , 










10.2 




> 


10 


2S ^\ 




* IP Z.ZZ 




9.8 






9.6 


is N 




9.4 






9.2 





0.2 0.4 0.6 0.8 1 

T/T c 

Figure 2. The masses of different quarkonia states and the open charm (beauty) thresh- 
old as function of the temperature. Shown are the charmonia masses and open charm 
threshold (top) and bottomonia masses and open beauty threshold (bottom) as function 
of the temperature. The thick solid line is the open charm (beauty) threshold obtained 
from normalization at r = 1/(4T). The thin solid line in the open charm (beauty) 
threshold obtained from normalization at r = \/2/(4T) (see text). 



V string (r,T) = -(e - iarctan(2rT)) 1 + 

7rT 2 2T 2 1 1 

(a(T) - — _ ar ctan— ) r + - ln(l + 4r 2 T 2 ). (5) 

This form of the potential describes quite well the temperature dependence of the heavy 
quark potential in quenched QCD for appropriately chosen <r(T) [17]. In order to make 
contact to real QCD we set e = 0.4 and use T c /^/a = 0.425 from [12] for the decon- 
finement temperature ( by a we always denote the string tension at zero temperature). 
Furthermore we use the following values of the heavy quark masses, m c = 1.3GeV and 
nit = 4:.72GeV as well as ifa = OAiGeV for the zero temperature string tension. This 
set of parameters gives a fairly good description of the observed quarkonium spectrum at 
zero temperature. The temperature dependence of the string tension was taken from [17]. 
The resulting quarkonia masses are shown in Fig 2. 

Since the smallest distance available on lattice is only 0.25T -1 one may worry about 
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possible medium effects at this distance and their role in determination of V^o(T). Nor- 
malizing the Polyakov loop correlator (2) at r — 1/ (4T) with Eq. (5), we thus obtain what 
might be a more reliable estimate of the plateau Voo(T) than with the T = form (3). It 
turns out, however, that the two forms of short distance behavior resulting from the zero 
temperature Cornell potential (3) and (5) are practically identical, so that the normaliza- 
tion is in fact not affected by the in-medium modifications at larger distances. To consider 
further possible uncertainties of the normalization procedure, the Polyakov loop correla- 
tor has been normalized also at the next smallest distance r = y/2/(4T). The resulting 
two forms of Vqo{T) are shown in both Fig. 2. The difference between the two curves of 
Voo (T) provides an estimate of the normalization error. Except for the region very near 
T = T c , the uncertainty is seen to be quite small. 

From Fig. 2 one can see that ip' and \c states become an open charm states well below 
T c and can dissociate by decaying into DD. The situation is similar for T(3S) and Xb(2P) 
states which can decay into BB below T c . For J/tp, Xb(lP) and T(25) it is not possible 
to say whether they will dissociate above T c or just below T c . Finally, the T(15) state will 
definitely dissociate above the deconfinement. We are going to estimate the dissociation 
temperatures of these states in the next section. 



4. Quarkonium dissociation by color screening and the sequential suppression pat- 
tern 

In the deconfined phase it is customary to choose the constant C in (2) to be the value of 
the correlator at infinite separation C = Tin < L(r)L t (0) > = In | < L > | 2 . The 

resulting connected correlator defines the so-called color averaged potential [18] 

v ( ,,D = -rh <^f"'»> (6) 

| < L > | 

The color averaged potential can be written as the thermal average of the potentials in color 
singlet V\ (T, r) and color octet Vg (T, r) states: 

exp(-V(r,T)/T) = ±e X p{-V 1 (r,T)/T) + | cxp(-y 8 (r, T)/T) (7) 

In potential models it is assumed that quarkonium is dominantly a singlet QQ state. Fur- 
thermore the octet channel is repulsive (at least in perturbation theory) and therefore only 
a singlet QQ pair can be bound in the deconfined phase. Thus we need to know the singlet 
potential. Lattice data in the relevant case of 3 flavor QCD exist only for the averaged 
potential [12,19]. The averaged potential in 3 flavor QCD is shown in Fig. 3 for three rep- 
resentative temperatures. Note that within the present accuracy of the lattice calculations 
the potential vanishes beyond some distance r {T) denoted by vertical arrows in Fig. 3. In 
perturbation theory, the leading terms for both singlet and octet potentials are of Coulomb 
form at high temperature and small r (r << T -1 ) 

V l{ r,T) = - 4 -^l, V s (r,T) = + 1 -^l, (8) 
6 r or 

with a{T) for the temperature-dependent running coupling. In the region just above the 
deconfinement point T = T c , there will certainly be significant non-perturbative effects 
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Figure 3. The color average potential calculated on lattice 

of unknown form. We attempt to parameterize the existing non-perturbative effects by the 
following form 

Vi(r,T) = -l^lexp{- M (T)r}, V s (r,T) = I C (T)^ exp{- M (T)r} 
6 r or 

(9) 

where /x(T) denotes the effective screening mass in the deconfined medium. We fit the 
lattice data of [12, 19] by Eqs. (7) and (9) assuming a(T) and )i(T) and c(T) to be unknown 
functions of T. Since the present data are not precise enough to determine a{T) and n{T) 
and c(T) simultaneously, some additional constraints coming from simulations of pure 
gauge theory [20-23] should be invoked in the fit procedure. The fit procedure is described 
in detail in Ref. [11]. Here it is sufficient to mention that a(T) can be well described by 
1-loop formula for the QCD running coupling with Aqcd — (0.34 ± 0.01)T C and the 
screening mass n(T) is constant in units of the temperature ^{T) = (1.15 ± 0.02)T [11]. 

Now we are in a position to discuss quarkonium dissociation due to color screening. 
It is natural to assume that the heavy QQ pair cannot exist as a bound state if its effec- 
tive binding radius ( the mean distance between Q and Q) is larger than the screening 
radius of the medium. The effective radii for different bound states are calculated from Eq. 
(4) with V(r) = Vi(r, T). The screening radius of the medium can be identified with 
l//i(T). However, the value of /i(T) strongly depends on assumptions we have made to 
determine it. A less model dependent and more conservative approach would be to identify 
the screening radius with r (T) defined above. We use the latter approach. In Fig. 4 we 
show the effective radius of J/ tp and T states and the screening radius as function of the 
temperature. The intersection of these curves defines the dissociation temperature of J/ tp 
and T states. Similar analysis was done for excited states which may survive above T c . 

Alternatively one can define the dissociation temperatures as the temperature where the 
effective bound state radius diverges [24] and the quark-antiquark pair is unbound. It is 
clear from Fig. 4 such definition will lead to larger, though not very different value of 
the dissociation temperature. However, one should note that once the radius of the bound 
state is larger than the screening radius the present treatment based Schrodinger equation 
is clearly not valid and the effect of the medium becomes so strong that quarkonium disso- 
ciation is very likely to happen. 
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Now the dissociation temperatures are known for all quarkonium states and summarized 
in Tab. 2. Combining these dissociation temperature with the feed-down fractions deter- 
mined in section 2 we can predict the sequential suppression pattern of J / ip and T states 
as function of the temperature. These are summarized in Fig. 5. 



qq 


J/* 


Xc 


V>' 


T(15) 


Xb(lP) 


T(25) 


Xb(2P) 


T(35) 


T/T c 


1.10 


0.74 


0.2 


2.31 


1.13 


1.10 


0.83 


0.75 



Table 2: Dissociation temperatures of different quarkonium states. 



5. Lattice determination of quarkonium properties at finite temperature 

At zero temperature quarkonium properties (masses , decay constants etc.) can be obtained 
from the behavior of the mesonic correlators at large Euclidean time separations 

<O(T,p),O(0,-p) >™=^A„e- m " T (10) 

n 
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here 0(r,p) = J^s e lpx O(r, x) is the mesonic operator bilinear in quark fields 0(t, x) — 
q(x)Tq(x); T = 75, 7^ for the pseudo-scalar and the vector channel, respectively. Direct 
application of Eq. (10) for determination of the quarkonium masses is difficult because 
at large separations statistical errors in the correlation function become very large. To 
improve the situation one usually uses smeared (improved) operators O instead of O such 
that only one coefficient A n in Eq. (10) is significantly different from zero. A possible 
choice of the operator O which has an optimal projection onto a given quarkonium state is 

6(r, z) = ^u(j>-(t, z)Tq(r, z + y), (11) 
y 

where ui(y) is the trial wave function 2 (see [25] for further details). 

At finite temperature further difficulties appear in direct lattice determination of quarko- 
nium masses. First of all Eq. (10) is not applicable in principle because the time extent 
is limited by the inverse temperature. Since the temporal lattice size becomes smaller as 
the temperature is increased fewer data points on the temporal correlators are available. 
Another problem which is also present at zero temperature but becomes more serious at fi- 
nite temperature (just because one is enforced to consider mesonic correlators also at short 
distances) is the discretization errors of order niQa with ttiq being the heavy quark mass 
and a being the lattice spacing. 

To study quarkonium properties at finite temperature in quenched QCD correlators of 
smeared operators were calculated on anisotropic lattices [26,27]. To deal with large dis- 
cretization errors in Ref. [26] the NRQCD formalism was used (where the scale mq is 
integrated out). In this calculation the mesonic correlators corresponding to the ground 
state quarkonium show only small changes up to temperature l.2T c compared to the zero 
temperature case while the first excited state shows very dramatic change with the temper- 
ature. The authors of Ref. [27] used a different formalism, namely they used the Fermilab 
action [28] on anisotropic lattices which has no discretization error of order rriQa at tree 
level. They observed more dramatic changes of the mesonic correlators as the deconfine- 
ment point is crossed, however, the correlators are very far from the free ones even at 
temperatures as high as 1.5T C indicating possible existence of bound states at this temper- 
ature. 

There are several problems with the approach presented in Ref. [26,27]. The most seri- 
ous one is the use of optimized correlators which assumes existence of well defined quarko- 
nium states at finite temperature. Clearly, different particle states at zero temperature will 
appear as quasiparticles with finite width at non-zero temperature. The implementation 
of NRQCD formalism in [26] does not assume anti-periodic boundary condition in time 
direction for the quark propagators and therefore it is not clear to what extent the meson 
correlators calculated in [26] can be related to some finite temperature (retarded) corre- 
lation functions. In Ref. [27] the spatial lattice spacing was quite large and the effect of 
doublers cannot be completely neglected. 

Another possibility to extract quarkonium properties in lattice QCD is to calculate 
usual point-to-point correlators G(r,p) =< O(r,p)O(0,—p) > (i.e. correlators of point 



2 This definition works only in Coulomb gauge 
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sources) in imaginary time and extract the spectral function. The information on the bound 
state are then encoded in the peaks of the spectral function. 

The imaginary time correlator can be related to the retarded meson propagator Gr(u,p) 
by analytic continuation 



In principle this equation allows to determine a(ui) from G(r,p). However, in practice 
the value of G is available only for N T /2 <~ 8 different r -values. For reasonably fine 
discretization in w -space one has <~ 700 degrees of freedom to be reconstructed. This 
problem can be solved only using the Maximum Entropy Method (MEM) (see [29] for 
review). Mesonic spectral functions in different channels were successfully reconstructed 
using this method [30]. Very recently it has been demonstrated that spectral function also 
at finite temperature can be reconstructed using this method and isotropic lattices with 
N T = 12, 16 sites in temporal direction [31,32]. Thus a natural alternative to the approach 
presented in Ref. [26,27] could be the calculation of point-to-point mesonic correlators 
using isotropic lattice with non-perturbatively improved clover action [32]. Compared 
to Ref. [27] this approach has the advantage that all lattice artifacts of order Ta pa are 
completely removed but the also the disadvantage that mqo effects are present even at the 
tree level. The latter, however, can be controlled if a is small enough. This approach is 
currently being investigated [33]. Calculations so far have been performed on 48 3 x 16, 
48 3 x 12 and 64 3 x 16 lattices and gauge coupling j3 = 6.499, 6.640 and 6.872. These 
parameters correspond to temperatures T = 0.9T C , 1.2T C and 1.5T C and lattice spacing 
from a -1 = 6.5GeV, 4.9GeV and 4.0GeV. The quark mass in these calculations was 
about l.hGeV. The discretization errors in the correlation function turns out to be large, e.g 
the value of G(t = 1/(2T)) at 1.5T C calculated with a" 1 = 4.9GeVanda~ 1 = 6.5GeV 
differs by 16% while the same difference is the light quark sector is about 1%. Because 
of this MEM cannot be applied to extract the spectral function unless extrapolation to the 
continuum limit is performed. Nevertheless, correlation functions themselves can provide 
some information about existence of quarkonium bound states. In Fig. 6 the ratio of the 
meson correlator calculated at 1.2T C (48 3 x 12, (3 = 6.499) to the one calculated at 0.9T C 
(48 3 x 16, j3 = 6.499). Except the last point at t ~ 0.3/m this ratio stays very close 
to unity at la level. The deviation at r ~ 0.3/m is simply due to periodic boundary 
condition (r ~ 0.3 fm corresponds to r = 1/(2T) on 48 3 x 12 lattice). While this 
fact does not necessary implies that the quarkonium spectral function does not change as 
the deconfinement point is crossed, it does imply that the propagation of heavy quarks in 
the deconfined phase is far from free propagation of QQ pair in agreement with studies 
performed in Refs. [26,27]. 
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which allows to write down the spectral representation for G(t, p) 
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Figure 6. The ratio of the correlators in pseudo-scalar and vector channels at 1.2T C to 
the one at 0.9T C . Calculations were done on at j3 = 6.499 on 48 3 x 12 and 48 3 x 16 
correspondingly. 



6. Conclusions 



I have considered quarkonium dissociation in hot strongly interacting matter below as well 
as above the deconfinement. In a confined medium dissociation of certain quarkonium 
states occurs due to in-medium modification of the open charm (beauty) threshold as well 
as the quarkonia masses. In the deconfined medium quarkonium dissociation is due to 
color screening. In this analysis quarkonium masses were extracted from Schrodinger's 
equation. The singlet potential used in Schrodinger's equation was extracted from lattice 
data on the Polyakov loop correlator using some additional assumptions. For a more accu- 
rate determination of the quarkonium suppression patterns, it would be desirable to carry 
out direct lattice studies of the color singlet potential and of its quark mass dependence, 
which may become important near the critical temperature. Furthermore, to make contact 
with nuclear collision experiments, a more precise determination of the energy density via 
lattice simulations is clearly needed, as is a clarification of the role of a finite baryochemi- 
cal potential. For the latter problem, lattice studies are so far very difficult; nevertheless, a 
recent new approach [34] could make such studies feasible. 

In section 5 the problem of direct lattice determination of quarkonium properties was 
discussed. An unbiased determination of these properties free of lattice artifacts seems to 
be very difficult and not available so far even in the quenched approximation. Such an 
analysis will be a very important check of the results obtained from the simple potential 
picture discussed in sections 3 and 4. 
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